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Roots of a Polynomial (SADOI Only) 
Closed subroutine 
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This subroutine calculates the roots of the polynomial 

n-2 . 



a * n + a l * n " 1 + °2 X 



The coefficients a Q , a. 



+ a n-l x + a n " ° 



. . . a are to be placed in 
' n 



*0* 1, 
the Williams Memory beginning at (COEFF). Enter the 



subroutine with Q = 50 ■ n 50 q. where n is the degree of 
the polynomial. This will cause complex roots to be 
found in conjugate pairs with a saving in timeo The 
entry Q = JO n 50 q. will cause the roots to be found 
individually. The latter entry should be used when 
the coefficients are not all real. 

STORAGE AMD SYMBOLIC ADDRESSES: 

The n + 1 coefficients should be placed with a Q in 
(COEFF), a x in 2(C0EFF), a 2 in MCOEFF), etc. Thus, 
2n + 2 words are necessary at (COEFF). 
The routine may put an extraneous root in 2n(R00TS) 
and 2n + 1 (ROOTS). Thus 2n + 2 words of storage are 
necessary at (ROOTS). 

2k words of temporary storage are required. This 
block is to be specified by S^. 
(A5) must define the complex number floating point 
arithmetic routine. See this routine for number ^ 
format, etc. 

S3 is the location of the floating accumulator used 
with A5. 
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IHffiATIOH: 



ACCURACYs 



MATHEMATICAL METHOD; 



B-UNES USED: 



0, 1, 2 S and 3 are used by A5 while this routine is 

operating. 

The time is approximately 3 seconds per root. This 

is reduced if complex roots are found in conjugate 

pairs. The time also depends on the condition of 

the polynomial. In any event, the time will he no 

greater than that of Library Routine J 2. 

Usually about 9 decimal places. Multiple roots are 

not found as accurately. 

A root R is found using the same method as J 2. 

yjPj ~ -ft 

If - U-p ' ^ l > ~ 10* the root is classified as complex. 

For the 50 n 50 q entry, the conjugate of R. is placed 
in the block of roots and divided outo For the 
JO n 50 q entry, the procedure is exactly that of J 2. 
B-lines and 1 are used during the operation of this 
routine and are not reset. 
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